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Abstract 

Theoretical understanding of the observed jet quenching measurements at RHIC and LHC is 
challenging in QCD because it requires understanding of parton to hadron fragmentation function 
in non-equilibrium QCD. In this paper, by using closed-time path integral formalism, we derive 
the gauge invariant definition of the gluon to hadron fragmentation function in non-equilibrium 
QCD which is consistent with factorization theorem in non-equilibrium QCD from first principles. 
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I. INTRODUCTION 


During the early stage of the evolution of the universe, just after ~ 10“^^ seconds of the 
big-bang, the universe was hlled with a state of matter known as quark-gluon plasma. This 
early stage of the universe is known as the quark epoch. The temperature of the quark-gluon 
plasma is very high ^ 3.2 x 10^^ K 200 MeV) which is about million times larger than 
the temperature of the sun (~ 1.56 x 10^ K). Besides black holes, the quark-gluon plasma 
is much denser than all the forms of the matter we know so far such as neutron stars etc.. 
Hence it will be a huge step for us if we can recreate this early universe scenario in the 
laboratory, i. e., if we can produce quark-gluon plasma in the laboratory. At present, the 
RHIC and LHC heavy-ion colliders are the best facilities to create quark-gluon plasma in 
the laboratory. 

RHIC (relativistic heavy ion collider) at BNL collides two gold nuclei at = 200 GeV 
and the LHC (large hadron collider) at CERN (in its first run) collides two lead nuclei at y/s 
= 2.76 TeV. Since the total energy at RHIC and LHC are very high ~ 40 TeV and ~ 574 
TeV respectively which are deposited in small volume at the initial time, the initial energy 
density might have been mnch higher than the eqnivalent temperature of ~ 200 MeV to 
produce quark-gluon plasma at RHIC and LHC. In the second rnn at the LHC in the PbPb 
collisions at = 5.5 TeV the total energy ~ 1150 TeV will even create much higher initial 
energy density. 

However, the main problem at RHIC and LHC is whether the quark-gluon plasma is 
thermalized or not. In case of a non-thermalized quark-glnon plasma at RHIC and LHC 
the notion of temperature does not exist and hence the hnite temperatnre lattice QCD 
stndies become inapplicable at RHIC and LHC. The qnestion of non-eqnilibrinm qnark- 
gluon plasma arises because the two nuclei at RHIC and LHC travel almost at speed of 
light. This means the partons inside the two nnclei just before the collision at RHIC and 
LHC carry very high longitndinal momentnm and very small transverse momentnm creating 
momentum anisotropy after the nuclear collision. In order for these anisotropic partons to 
form thermalized qnark-glnon plasma many more secondary partonic collisions over a large 
period of time is required. However, the typical hadronization time in QCD is very small 
(~ seconds) after which the partons become hadrons. Hence there may not be enough 
time for these anisotropic partons to thermalize before hadronization takes place at RHIC 
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and LHC. Hence it is necessary to stndy noneqnilibrinm-nonpertnrbative QCD at RHIC and 
LHC which, of conrse, is a difhcnlt problem. 

It shonld be remembered that the qnark-glnon plasma can not be directly detected at 
RHIC and LHC becanse qnarks and glnons are not directly observed. The primary indirect 


dnon plasma detection are, 1) j/^snppression |l|, 2) dilepton and direct 
2|, 3) strangeness enhancement |3| and, 4) jet qnenching etc.. We 


signatnres of qnark 
photon prodnction 
will focns on the jet fragmentation in non-eqnilibrinm QCD in this paper. 


Non-eqnilibrinm qnantnm field theory can be stndied by nsing closed-time path (CTP) 
integral formalism ja ITJ. Hence in order to stndy noneqnilibrinm-nonpertnrbative QCD 
it is necessary to nse path integral formnlation of non-eqnilibrinm QCD by nsing CTP 
formalism. It shonld be mentioned, however, that implementing CTP in non-eqnilibrinm 
at RHIC and LHC is a very difhcnlt problem, especially dne to the presence of glnons in 
non-eqnilibrinm and hadronization etc. Even the calcnlation of one-loop glnon self energy 
in non-eqnilibrinm in covariant gange which is derived in 8|, |9|] becomes mnch more tedious 
than the calculation of one loop QCD diagram in vacuum due to the presence of additional 
closed time path indices in non-equilibrium. 

Recently, by using Schwinger-Keldysh closed-time path integral formalism, we have de¬ 
rived the gauge invariant dehnition of the quark to hadron fragmentation function in non- 
equilibrium QCD [lO, 11|. In this paper we will extend this study to derive the gauge 
invariant dehnition of the gluon to hadron fragmentation function in non-equilibrium QCD 
which is obtained from the single gluon m-state \g > in non-equilibrium QCD. 

The paper is organized as follows. In section H we present the main results and discuss 
the novel features. In section HI we give a brief description of closed-time path integral 
formalism in non-equilibrium QCD. We implement closed-time path integral formalism in 
non-equilibrium QCD to derive the dehnition of the gauge non-invariant gluon to hadron 
fragmentation function in non-equilibrium QCD in section IV. In section V we prove fac¬ 
torization of infrared divergences in non-equilibrium QCD. In section VI we derive the 
dehnition of the gauge invariant gluon to hadron fragmentation function in non-equilibrium 
QCD. Section VH contains conclusions. 
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II. MAIN RESULTS AND NOVEL FEATURES 

Before presenting the technical details of the derivation of the gluon to hadron fragmenta¬ 
tion function in non-equilibrium QCD we will hrst summarize the main result and its novel 
features in this section. The readers who are not interested in the technical details can skip 
the rest of the paper. 


A. Main Results 


For the gluon with arbitrary non-equilibrium (non-isotropic) distribution function fg{k) 
at initial time to, we hnd that the gauge invariant dehnition of the gluon to hadron frag¬ 
mentation function in non-equilibrium QCD which is obtained from the single gluon m-state 
\g > and is consistent with factorization theorem in non-equilibrium QCD is given by 




jd-2 


Dh/s{z,Pt) ife + ( 2 ^)d-i 

< in\Q^'"{x~, xt)^[x~,XT] a^HiP~^^ 0T)aH{P^, 0r)<h[0]Q 


^ik'^x -\-iPT ■Xt/z 


in > 


(1) 


which is valid in covariant gauge, in light-cone gauge, in general axial gauges, in general 
non-covariant gauges and in general Coulomb gauge etc. where Q^°‘{x) is the (quantum) 
gluon held which fragments to hadron and |m > is the initial state of the non-equilibrium 
QCD medium in the Schwinger-Keldysh in — in closed-time path formalism and the path 
ordered exponential 


^ab\x\ = 


POO 

V exp[-7^ / d\l- A"(x + /A) 


rp{-A)c _ . rcab 

^ab — 


( 2 ) 


*^0 J ab 

is the non-abelian phase or non-abelian gauge link in the adjoint representation of SU(3) 
where is the light-like four-velocity and is the SU(3) pure gauge background held 

with a,h,c = 1,2,..., 8. The dehnition of the gluon to hadron fragmentation function in 
non-equilibrium QCD in eq. ([T]) is gauge invariant with respect to the gauge transformation 


T^A'Ax) = U{x)rAl{x)U-\x) + - [dgU{x)]U-\x) 


W 


along with the homogeneous transformation 12l4l5| 


T^Q'^{x) = U{x)T^Q‘;,{x)U-\x) 


( 3 ) 

( 4 ) 
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where 


U{x) = 


(5) 


The special case fg{k) = — corresponds to the hnite temperature QCD in equilibrium. 

e^-l 

We will present a derivation of eq. ([1]) in this paper. 


B. Novel and Qualitatively Interesting Features 


Let us describe the novel features of the dehnition of the gluon to hadron fragmentation 
function in non-equilibrium QCD obtained in eq. O- 

Similar to DGLAP evolution equation 1^ in pp collisions we will follow the DGLAP-like 
evolution equation procedure in AA collisions at RHIC and LHC where the inclusive hadron 
production cross section is given by the factorized formula in non-equilibrium QCD 


dydp^ 


_ f dz (i(Jj 

\ J dydpli 




( 6 ) 


where is the partonic level differential cross section in non-equilibrium QCD and 

Du/ii^iQ) is the fragmentation function in non-equilibrium QCD. The evolution of the 
fragmentation function DHpiz, Q) in non-equilibrium QCD obeys the generalized DGLAP- 
like evolution equation in non-equilibrium QCD 


dlnQ Ti ^ Jz z' z' 


(7) 


where Pi^jk{z',Q) is the splitting function in non-equilibrium QCD 17|]. In the above 
equations the indices i,j, k = q,q, g. 

Note that one of the crucial ingredient to derive DGLAP evolution equation is the factor¬ 


ization theorem 


16| . Hence similar to DGLAP equation in QCD in vacuum the DGLAP-like 


equation in non-equilibrium QCD in eq. ([7]) is consistent with the factorization theorem in 
non-equilibrium QCD (see also eq. (jb]) and appendix A). This implies that the partonic 
level differential cross section in eq. ([b]) in non-equilibrium QCD and the fragmenta¬ 

tion function DH/i{z,Q) in eq. ([6]) in non-equilibrium QCD are studied by using \in > for 
the ground state in non-equilibrium QCD instead of the usual vacuum state |0 > in QCD. 
It should be remembered that the ground state \in > in non-equilibrium QCD contains 
both the vacuum part and the medium part. For example, when the distribution function 
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f{k) = 0 then \in > becomes |0 > and one reproduces all the equations and quantities of 
the QCD in vacuum. 

The leading order perturbative gluon propagator in non-equilibrium QCD is given by 




= -ils"- + (a - l)-p-] GT:{k) - 


( 8 ) 


which contains the vacuum propagator plus the medium propagator where 


and 


tAt. -2^s(e)e(- 

[ -27,s(k^)e(k„) -pi- 

. . -CD 

G“y(^ = 2irS(k:‘)f,(k) 



(9) 


( 10 ) 


In the above equations the closed-time indices r, s = -f, — correspond to upper and lower 
time branch in the closed-time path formalism. 

However, since the gluon fragmentation function is a non-perturbative quantity in QCD 
it is not possible to decompose the gluon fragmentation function into the vacuum part and 
the medium part (see eq. ([T])). Hence one finds that to be consistent with factorization 
theorem in non-equilibrium QCD, the fragmentation function in QCD in vacuum can not 
be used in the DGLAP-like equation in non-equilibrium QCD in eq. ([7]) to study hadron 
production from quark-gluon plasma. The fragmentation function in non-equilibrium QCD 
should be used in DGLAP-like equation in non-equilibrium QGD in eq. ([7]) to study hadron 
production from quark-gluon plasma by using factorized formula as given by eq. (E]) where 
is the partonic level differential cross section in non-equilibrium QGD and Dn/ii^z, Q) 
is the fragmentation function in non-equilibrium QGD. 


III. NON-EQUILIBRIUM QCD USING CLOSED-TIME PATH INTEGRAL FOR¬ 
MALISM 

When a system is in non-equilibrium its asymptotic future is quite different from its 
initial preparation in the remote past. Hence unlike zero-temperature field theory or finite 
temperature field theory in equilibrium there is no simple relation between asymptotic future 
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state and initial state in the remote past. For this reason the generic non-eqnilibrinm 
qnantnm held theory is very non-triviah In order to avoid the asymptotic fntnre state and 
to deal with only initial state in the remote past one introdnces two time branches in the 
closed-time path (CTP) formalism 6|, l7|. 

The glnon to hadron fragmentation fnnction in non-eqnilibrinm QCD in eq. ([T]) is a 
non-pertnrbative qnantity. Hence it is necessary to nse path integral formnlation of non- 
eqnilibrinm QCD in CTP formalism to stndy its properties. Before going to glnons in 
non-eqnilibrinm QCD let ns consider the scalar glnons hrst in the CTP formalism in non- 
eqnilibrinm in the path integral formnlation. The generating fnnctional in scalar held theory 
in non-eqnilibrinm is given by 


Z\p, J+, J_] = J exp|i J d*x[L[(l>^] - 0|/3|0, > 


( 11 ) 


where p is the initial density of states in non-eqnilibrinm, £[0] is the fnll lagrangian density in 
scaler held theory and |0±, 0 > is the qnantnm state corresponding to the held conhgnration 
= 0) at the initial time. The snbscript -|- refers to closed-time path index in the 
positive time contonr and the snbscript — refers to closed-time path index in the negative 
time contonr in the closed-time path formalism 6|, l7|. 

Since there are two time branches — in the closed-time path formalism there are two 
external sonrces J+, J_ corresponding to these two diherent time branches. Hence one hnds 
that there are fonr green’s fnnctions in non-eqnilibrinm which are given by 


G++{x,x') = =< T(t){x)(t){x') >=< m|T0(a;)0(a;')|m > 

i^oJ+[x)J+[x') 

=< >=< tn\Td>(x)^(x’)\m > 

[—iydJ_[x)J_[x') 

G+-{x,x') = w’ (t){x')(t){x) >=< in\(t){x')(t){x)\in > 

G-+{x,x') = tAw’ (t){x)(t){x') >=< in\(t){x)(t){x')\in > (12) 

where \in > is the initial state in non-eqnilibrinm qnantnm held theory at the initial time 
t = tin = 0. Note that dne to the presence of the medinm the initial state \in > in non- 
eqnilibrinm qnantnm held theory is diherent from the vacnnm state |0 > in the qnantnm 
held theory in vacnnm. 
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In order to extend scalar glnon to gluon in non-equilibrium QCD we need to add gauge 
fixing term in the lagrangian density which results in the appearance of the Faddeev-Popov 
(F-P) determinant in the generating functional which can be expressed in terms of path 


integral over the ghost helds 


18| . However, we will directly work with the Faddeev- 


^opov 
91 for 


(F-P) determinant in this paper and we will work in the frozen ghost formalism 
the medium part at the initial time t = tin = Hence it is now straightforward to extend 
the generating functional of the scalar gluon in non-equilibrium in eq. flTT]) to gluon in 
non-equilibrium QCD. 

Extending the generating functional of the scalar gluon in non-equilibrium in eq. flTT]) 
to gluon in non-equilibrium QCD we find that the generating functional in non-equilibrium 
QCD in the path integral formulation is given by |8|, l9| 


Z[p, J+, J_, T]+, r]_, r/+, r/_] = 
1 




exp 


] [dQ-] [#+] [#-] [#+] [#-] det( )det( ) 

[i I - F‘liQ.]) - 

+ J+ ■ Q+ - J- ■ Q- + f]+ ■ ^l^+ - v- ■ + ^p+ ■ ri+ - ■ ri_]] 

X < Q+,'0+,^+,O| P > (13) 


where p is the initial density of state in non-equilibrium QCD, the state |Q^, "0^, 0 > 

corresponds to the field configurations Q))(F, t = tin = 0), = tin = 0), t = Un = 0), 

the J^^{x) is the external source for the quantum gluon held Q^'^{x), the f]i{x) is the external 
source for the Dirac held ^pi{x) of the quark, a is the gauge hxing parameter, det(^^^^^±^) 
is the Faddeev-Popov (F-P) determinant and 


F^Q±] = d.Qt^ix) - a.Q^±(x) + pr“Qt±(x)Q“±(x), 

I/[Q±] = *7'^5;. + pTVQm±- (14) 


(5 (d ) 

As mentioned above the Faddeev-Popov (F-P) determinant det(—£y^) in eq. flT^ can be 


expressed in terms of path integral over the ghost helds m but we will directly work with 


the Faddeev-Popov (F-P) determinant det( 




5a;S 


) in this paper. 


Note that the generating functional in non-equilibrium QCD in eq. flT3l) is generic from 
which the hnite temperature QCD formalism and QCD in vacuum formalism can be ob¬ 
tained. 

From eq. flT^ we hud that the non-perturbative gluon correlation function in non- 














equilibrium QCD is given by 


< in\Q^^rixi)Qlsix2)\in >= J [dQ+][dQ-][d^l)+][d^l)-][d^l)+][d^l)-] Q1 ^^{xi)QI^{x2) 


dd 

det CTl^ )det( )exp[i 




[Qh 


5ui ' ' 5ui 4^ ' 4^ 2a 

+'0+^[Q+]'0+ -^-^[Q-]'0-]] < Q+,'0+,^+,o| p |o,'0_,'0-,Q- > 




(15) 


where r, s = +, — are the closed-time path indices (” -f ” sign corresponds to positive time 
oath and ” — ” sign corresponds to negative time path in the closed-time path formalism 




IV. GAUGE NON-INVARIANT DEFINITION OF GLUON TO HADRON FRAG¬ 
MENTATION FUNCTION IN NON-EQUILIBRIUM QCD 

For simplicity, let us consider the scalar gluon in non-equilibrium hrst before considering 
the gluon in non-equilibrium QCD. For the application to RHIC and LHC high energy 
nuclear collisions it is useful to consider the light-cone coordinate system x^ = {x'^,x~,xt)- 
Note that at the initial time x'^ = xf^ = 0, the interaction picture coincides with the 
Heisenberg picture and Schrodinger picture. Hence at the initial time x~^ = xf^^ = 0 we 
write the scalar held as 

( 16 ) 

where a^{p) and a{p) are the creation and the annihilation operators of the scalar gluon 
respectively. The commutation relations at the initial time are given by 

Hp^,PT),a\p'^,PT)]x+=o = (2vr)'^"^h(p+ -Pt)i 

Hp^,PT),a{p^,PT)]x+=o = WiP'^,PT),a\p'^,PT)]x+=o = 0. (17) 

The non-equilibrium distribution function f{p'^,PT) of the fragmenting scalar gluon at initial 
time is given by 

< in\a^{p'^,pT)a{p'^,PT)\in >= /(p+,pr)(2vr)'^“^5(p+ - p'^)S^‘^~‘^\pt - Pt) (18) 

where we have assumed the space {x~ and xt) translational invariance at initial time x~^ = 
xf^ = 0. The special case fg{p) = — corresponds to the hnite temperature held theory 

in equilibrium. 
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At the initial time x'^ = xf^ = 0 the scalar gluon state in the non-equilibrium quantum 
held theory is given by 

\p'^,Pt >= a\p'^,pT)\in > (19) 

where a^{p) is the creation operator of the scalar gluon and \in > is the initial state of the 
non-equilibrium medium. By using eqs. flT^ . flT8|) and flT7|) we hnd, at initial time, 

< p~^,P t\p~^,P T >=< in\a{p~^,pT)a\p~^,p rp)\in > 

= (27r)‘^"^5(p+ -p't)[1 + fiP^^Pr)]- (20) 

Consider the inclusive production of hadron H created in the out—state \H + X > from 
a scalar gluon in non-equilibrium in the initial state \p > with the probability amplitude 

<H + X\p>, (21) 

where X being other outgoing hnal state hadrons. The distribution hp{P) of the hadron H 
with momentum P from the parton of momentum p in non-equilibrium can be found from 
the amplitude in eq. (1^ by using 

^ <P,Pt\H + X >< H + X\p-^',PT>=hp{P) <p-^,Pt\p'^\pt> ■ (22) 

We write the left hand side as 


^ <P+,Pt\H + X><H + X\p+',p'^>=J 2 <P^,PT\aUP)\X><XMP)\p^’,PT> 

X X 

=< p+,pT\a^HiP)aH{P)\p^', p't > • ( 23 ) 


From eqs. fl22|) . fl2^ . flT^ and fl20|) we hnd 

< in\a{p* ,pT)a}[i(P*<P t)iiu{P*,PtW{j>'*, p'T)\i'n > 

= 2z(2i,Y-^Du,,iz,PT) \l + f(p*,PTmp*- p'*)S‘‘-Ypt-p't)- ( 24 ) 


From eq. ffTOj) we obtain 


(27r)'^ ^<in\(t){x ,XT)a\{P^ ,PT)aH{P^ ,PT)(l){^)\in >= J ^Pt 

[< in\e~''P''^a{p^,pT)a'^H{P^,PT)aH{P^,PT)a^{p'^,p'T)\'^n >],e+=o- 




( 25 ) 
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Using this in eq. (|2^ we find the expression of the scalar gluon fragmentation function 
in non-equilibrium 


< m|0(x",a;r)a]:^(P+,Pr)aH(^^,^T)0(O)|m > . 


d 


d-2 


XT 


ip'^x —ipT'XT 


( 26 ) 


■ 3v making a Lorentz transformation to the zero transverse momentum frame of the hadron 
10l | we find from eq. fl26|) that the scalar gluon fragmentation function in non-equilibrium 


is given by 


Dniaiz, Pt) = 


p-1 


^ [1 + /(p+,Pt)] 


dx~ 


d 


d-2 


(27,) 


^ip+x ■^iP’j’-X'j’!z 


d-1 


< in\(f){x , a;r)a]^(P^, Or)a_H'(P^, O'r)0(O)|m > . 


( 27 ) 


In the above expression, f{p'^,PT) is the non-equilibrium distribution function of the frag¬ 
menting scalar gluon at initial time x^ = = 0 and \in > is the initial state of the 

non-equilibrium medium in the Schwinger-Keldysh in — in closed-time path formalism. 

In analogous to eq. fl^ for the scalar gluon case, we find the gauge non-invariant 
definition of the gluon to hadron fragmentation function in non-equilibrium QCD is given 
by 


Dh/,{z,Pt) ife [i + /^(p+_p^)]y * (2„)^-i 

< in\Q^''{x~,XT)a^HiP^7^T)aH{P^,0T)Ql{0)\in > 


d 


d-2 


^ip'^x -\-iPT-XTIz 


( 28 ) 


where /g(p^,pr) is the non-equilibrium distribution function of the fragmenting gluon in 
QCD at initial time, Q^^{x) is the (quantum) gluon field in non-equilibrium QCD which 
fragments to hadron H and \in > is the initial state of the non-equilibrium QCD medium 
in the Schwinger-Keldysh in — in closed-time path formalism. 

Note that the transition from massless scalar gluon fragmentation function in eq. (HT} 
to gauge non-invariant gluon fragmentation function in QCD in eq. fl28|l is achieved by 
the following procedure. The massless scalar gluon fragmentation function in quantum field 
theory in eq. fl27)) is proportional to the two-point scalar gluon correlation function of the 
type < in\((){x)a^fpaH(f){0)\in >. Note that, in quantum field theory, the two-point scalar 
gluon correlation function in scalar field theory is given by < in\(j){x)(j){0)\in > and the 
two-point gluon correlation function in QCD is given by < m|Q“(a;)Qy(0)|m > where 
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are Lorentz indices and a, b are color indices of the gluon fields Q^(a;) and Ql{0) in QCD. As 
mentioned above, in quantum field theory, the gluon fragmentation function is proportional 
to two-point gluon correlation function. Hence the gluon fragmentation function in QCD in 
eq. fl25D is obtained from the massless scalar gluon fragmentation function in eq. fl27|) by 
replacing 

< in\(j){x)a^jjaH<P{0)\in > ^ < in\Q'^{x)a^fjaHQ'^°‘{0)\in > ■ (29) 

Since the fragmentation function in QCD is a Lorentz scalar and color neutral quantity, the 
sum over Lorentz and color indices are taken in eq. fl28l) . The additional factor A eq. 
fl28|l is supplied because the average over the spin and color ( 2 ^) of tho gluon in QCD is 
taken which is absent in the scalar gluon case. 


V. PROOF OF FACTORIZATION OF INFRARED DIVERGENCES IN NON¬ 
EQUILIBRIUM QCD 


Since the gluon to hadron fragmentation function in non-equilibrium QCD in eq. fl28ll is 
a non-perturbative quantity we will use path integral formulation of quantum held theory 


for this purpose 0,1201. 

The generating functional in the background held me 


QCD in the path integral formulation is given by 


rod of QCD in non-equilibrium 


12l-ll5| 


Z[p, A, J+, J_,ri+,ri_,r]+,r]_] = J [dQ+][dQ_][dij+][d^_][d'il)+][d'ip_] 

X det( ^^^^*^~^ ) x exp[z J d^x[-^F^l^[Q+ + A+] + + A_] 

——(G"^(Q+))^ + + '^+^[Q+ + ^+]'0+ 

Za 2a 

+ A-]i)- + J+ ■ Q+ - J- ■ Q- + f]+ ■ ^l)+ - fj_ ■ ■ r]+ - ■ //_]] 

<Q+ + A+, -0+, -0+, 0| p |0, -0-, Q- + A_ > (30) 


where the gauge hxing term is given by 

G%Q±) = = ■D,.['4±]C± (31) 

which depends on the background held A^'^^x) and 

F',AA± + QA = dAAl± + Ql±] - dAAU + Q'±] + gr‘"lAl± + QU\AI± + «±] (32) 
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where is the quantum gluon held and is the background held 12l-ll4|. 

From eq. fl30|) we hnd that the nonequilibrium-nonperturbative gluon correlation function 
in the background held method of QCD is given by 


< 


in\Q^^^{xi)Ql,{x2)\in >a= j [dQ+][ciQ-][dt/’+][#_][#+][#_] QlrMQlsi.X2) 




^ )exp[i J d^x[--F‘^l^[Q+ + kl+] + -F^lAQ- + ^-] 

——(G'^(g+))^ + '^+^[Q+ + ^+]t/’+ 

+ A_]i)_]] < g+ + A+,-0+,-0+, o| p |o,0_,0_, g_ + > . 


For a light-like quark attached to inhnite number of gluons we hnd the eikonal factor [2 
d^k I ■ A'^ik) 2 f d^ki d^k2 TH ■ A^{ki)TH ■ A\k2) 


(33) 


i + grj 


le J 


+ ... 


(34) 


(27r)‘^ I ■ k + ie J {^tiY (2vr)^ {I ■ ki + ie){l ■ {ki -|- ^ 2 ) + 

poo 

= Pexp[f^/ dA/-kF(/A)Td 
Jo 

which describes the infrared divergences arising from the inhnite number of soft gluons 
exchange with the light-like quark of four-velocity F where V is the path ordering and the 
gluon held A^^{x) and its Fourier transform A^^{k) are related by 

r d^k 

A^Yx) = J 

The light-like quark traveling with light-like four-velocity F produces SU(^ pure gauge held 


Ak-x 


(35) 


A^'^{x) both in classical mechanics 2ll-l23| and in quantum held theory 2^ at all the time- 
space position x^ except at the position x perpendicular to the direction of motion of the 
quark (/ • x = 0) at the time of closest approach (xq = 0). When A^^^[x) = A^^{\1) as in eq. 
fl3T|) we hnd l-x = \l-l = \Y^ which implies that the light-like quark hnds the gluon held 
^Mc(x) in eq. fl3T|) as the SU(3) pure gauge. The SU(3) pure gauge is given by 


7’“3l'‘“(i) = -l|y4(i)|4>-‘(2^)l« 


(36) 


which gives the non-abelian gauge link [20 1 

poo 

$. .(a:) = [Pexp[-i^ / dXl ■ A%x + (37) 

0 

The gauge hxing term in eq. fl30l) [where G“(g) is given by eq. fl3l|l ] is invariant 

for gauge transformation of AF 


6Al = gr^'^Alu;^ + d,u‘^ 


(38) 
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provided one also performs a homogeneous transformation of HQ- 

SQi = gf‘^Q'‘y- (39) 

By changing the integration variable Q ^ Q — A in the right hand side of eq. fl33l) we hnd 
< in\Ql^{xi)Ql,{x2)\in >a= j [ciQ+][ciQ-][di/’+] [#_][#+][#-] 


^b \ Ab 


[Qt,{x2) - Al^{x2]detC-^^^)det{^^j^-^)exp[i J d^x[ 


■\KAQ+] + \K^[Q-] 


1-1 


^ + ^+^\.Q+^^+ - V’-^[Q-]^-]] < Q+, V'+, V’+, o| p |o,V’-, Q- > 


2a 


(40) 


where from eq. (l3Tll we hnd 

G%Q±) = a,QT + = D^IA±]Q^^ - d^A'S 

and from eq. (I39|) [by using eq. (138|) ] we hnd 

5Ql^ = -gr'^-u:lQl^ + d,u:l. 


(41) 


(42) 


Changing the integration variable from unprimed variable to primed variable we hnd from 
eq. dlO]) 

< i^\Qlr{^l)Qls{^2)\in >A= J [dQV] 

[Qt{x2) - A'1 ^{x2] X det(^^^^)det(^^^^)exp[i J 


2a 


(43) 


This is because a change of integration variable from unprimed variable to primed variable 
does not change the value of the integration. 

The equation 


Q"±i^) = Q^,±i^) + 17r“^±(^)Q,±(^) + d.uUx) 


(44) 


in eq. (142|) is valid for inhnitesimal transformation (a; << 1) which is obtained from the 
hnite equation 


T‘Q'‘Ax) = UAx)rQlj,{x)UZ\x) + -ld^UAx)]Ui\x), 

'^9 


u^{x) = 


( 45 ) 
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Simplifying infinite numbers of non-commuting terms we find 


C±(^) = + I 


^gM±{x) _ 

gM±{x) 


]cb [d^u:)^{x)]. 


where 


(x) = r^ojl{x). 


The fermion held transforms as 


(46) 


(47) 


i/.yi) = e«^'"l'"'>i/.±(i). (48) 

Under the finite transformation, using eqs. (I 4 fill . (14811 and we find 

Wil = W±l. l#ilWLl = W±]W±]. -PTlOil = = (a,Q^^(x))\ 

(49) 


i(4ln'‘9a - >H + 9 rVQ"±]t/’± = fcln'Se - m + 9 rVQ()±]t/’±. 


det[ 




= det[ 


Kd,Q^±{x)), 


Q':±((t) - y±((t) = ie'’“*‘'>]si,Qt±((t). 


(5a;± Suj± ’ ■ ■ 

Since we are working in the frozen ghost formalism at initial time the < 
Q+,'^+,'0+, 0| p |0,-0-, Q- > corresponding to initial density of states in eq. flT5|l is 
gauge invariant by dehnition. Hence from eqs. (H6|) and fHHj) we hnd [n| 

(50) 


< (5+,V’+,V'+,o| p >=< (5+,V’+,^+,o| p 10,-0-, > • 


Using eqs. and (1501) in eq. flT0]l we hnd 

< in\QlT{xi)QlXx2)\in >A= J [dQ+][dQ-][#+][#-][#+][#-] [e^^^^''^^]acQgr{xi) 

[ggMa(x 2 )]^^gd^(^^)det( )det( )exp[z J d*x[--F^l^[Q+] + 

^ )^ + V'+^[Q+]0+ < Q+,0+,0+,o| p |o, 0 _, 0 _,Q_ > . 

(51) 

Using the similar technique as above we hnd 

< >yi= J [dQ+][dQ_][d0+][#_][#+][d0_] 


AiO riiO (1^ r 

Qlrixi) Ql,{x2)det{—^^)det{—0—)exp[i J d^x[ 


doj^L 


?c2 
4 - gu 




^ + ' 0 +^[Q+] 0 + < Q+,'0+,0+,o| p |o,0_,0_,Q_ > 

(52) 
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where is given by eq. fHTD . From eqs. ffT^ and (1521) we find 


< in\Q‘l^^{xi)Ql^{x2)\in >=< >a (53) 


which is valid in covariant gauge, in light-cone gauge, in general axial gauges, in general non 


covariant gauges and in general Coulomb gauge etc. [20j. Note that the creation operator 
al and annihilation operator a of the quark are related to the quark field via the equation 
25| 


= EE 


spin p 


I m 

VEr 


[a{p)u{p)e -I- a){p)v{p)E'^'^ 


(54) 


where color indices are suppressed. Hence one finds that the quark field ■^(x) depends 
on the the creation (annihilation) operator (a) of the quark but is independent of the 
creation (annihilation) operator (o-b) of the hadron. Similarly the gluon field Q^°‘{x) is 
also independent of the creation (annihilation) operator {an) of the hadron. Since 
is independent of 'ip{x) and Q^°‘{x) one can perform exactly the similar steps of the path 
integral calculation as above to find from eq. fl53D that 




(55) 


which proves factorization of soft (infrared) divergences at all order in coupling constant in 
non-equilibrium QCD where [see eqs. fl47j) and fl37)) ] 


^bcr (•^) 


POO 

[Pexp[-z^ / dXl ■ H“(x + /A)T(^)“]]6„ (T(^)“)6, = 

Jo 


( 66 ) 


is the non-abelian gauge link or non-abelian phase in the adjoint representation of SU(3), E 


is the light-like four velocity and r = 


in the closed-time path formalism 


- corresponds to the upper and lower time branches 
l7|| and a, 6, c = 1, 2,..., 8 are color indices. 


VI. GAUGE INVARIANT DEFINITION OF GLUON TO HADRON FRAGMEN¬ 
TATION FUNCTION IN NON-EQUILIBRIUM QCD 

The definition of the gluon to hadron fragmentation function in non-equilibrium QCD in 
eq. (|28|) is not gauge invariant. To make it gauge invariant and consistent with factorization 
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of infrared divergences we need to incorporate Wilson lines. The non-abelian gauge link in 
the adjoint representation of SU(3) is given by the path ordered exponential, see eq. fl56|l . 

fOO 

= [Vexp[-tg / dXl ■ A%x + (57) 

Jo 

which under non-abelian gauge transformation, as given by eq. ([3]), transforms as 

= Uac{x)^cb{x) (58) 


where 

U^(x) = 


(59) 


and from eq. (jl]) we find 


Qtr(x) = 


(60) 


Hence by extendin g th e path integral formulation of quark fragmentation function in 


non-equilibrium QCD ll| to gluon fragmentation function in non-equilibrium QCD we hnd 
from eqs. fl55|) . fl56|) . fl28|) . flT^ . fl58|) and fl60l) that the gauge invariant dehnition of the gluon 


to hadron fragmentation function in non-equilibrium QCD which is obtained from the single 
gluon m-state \g > and is consistent with factorization theorem in non-equilibrium QCD is 
given by 


DH/g{z, Pt) — 




I62; [1 + fg{k+,kT)] 


dx~ 


d 


'd-2 


(27r 


^ik'^x +iPT ■XT IZ 


\d—l 


<in\Q^''{x ,xt)^[x ,XT]a];^(P'^,0T)af/(P“^,0T)$[0]Q^(0)|m > 


( 61 ) 


which reproduces eq. ([T]) which is valid in covariant gauge, in light-cone gauge, in general 
axial gauges, in general non-covariant gauges and in general Coulomb gauge etc. where the 
path ordered exponential is given by eq. ([2]). 

This completes the derivation of gluon to hadron fragmentation function in non- 
equilibrium QCD. 


VII. CONCLUSIONS 

Theoretical understanding of the observed jet quenching measurements at RHIC and 
LHC is challenging in QCD because it requires understanding of parton to hadron fragmen¬ 
tation function in non-equilibrium QCD. In this paper, by using closed-time path integral 
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formalism, we have derived the gauge invariant dehnition of the gluon to hadron fragmen¬ 
tation function in non-equilibrium QCD which is consistent with factorization theorem in 
non-equilibrium QCD from hrst principles. 
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